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ABSTRACT 


A three dimensional viscous marching analysis for supersonic inlets has 
-scent ly been developed. To verify this analysis several benchmark ax i sym- 
metric test configurations have been studied and are compared to experimental 
data. Detailed tno-dimensional results for shock-boundary layer interactions 
are presented for flows with and without bowdary layer bleed. T hree - 
dimensional calculations of a cone at angle of attack and a full inlet at 
angle of attack are also discussed and evaluated. Results of the calculations 
demonstrate tr» code's ability to predict caiplex flow fields and establishes 
guidelines for future calculations using similar codes. 


iNTROOUCTIOR 


Over the past few years, RASA Lewis ?e search Center has been increasingly 
involved in developing three-dimensional spatial marching analyses for the 
design of supersonic inlets. In support of this work, research was also 
initiated to systematically verify these computer analyses by a two tier veri- 
fication process. The first tier involves establishing the ability of the 
computer analyses to predict flow parameters of interest in "generic* con- 
figurations which would be typical of real world designs. A second tier veri- 
fication process was also initiated where the essential elements of the com- 
puter analyses are validated in very simple environments which highlight one 
or more basic cr "benchmark" flow mechanisms. This is called "benchmark" 
verification of the computer code. 

A thorough verification process for the PEPSIS computer code is being 
ccnd- cted at NASA Lewis. PEPSIS is a three-dimensional viscous marching 
ana ^sis for supersonic inlets of ax i symmetric, rectangular or conformal 
geometry which was jointly developed by Scientific Research Associates and 
NASA Lewis. Details of the algorithm and code development are given in 
(refs. 1 to 2). Some of the generic verification results have been presented 
in a paper by Anderson and Towne (ref. 3) and at the Lewis Computational Fluid 
Mechanics Workshop in 1980 (ref. 4). Anderson and Benson (re*. 5) present 
ome of the benchmark results for rectangular geometry verification, while 
this paper presents results for axisynmetric geometries. Four major test 
cases will be considered in this study; the Lewis 40-60 inlet (ref. 6) at zero 
angle of attack, the two-d imensior.il shock-boundary layer interaction of Rose 
(ref. 7), the cone at angle of attack by Rainbird (ref. 8), and the Lewis 
40-60 inlet at angle of att>_x. 
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AMLYSiS 


WSIS, * acr omy for Pa Aolk Elliptic StreMrise Implicit Svpersanic, 
is a fully tt re e d imensional copater de$ipi code for sapersonic inlets. The 
aaljsis retains the jeneral three-dimensional viscous nature of the lavier- 
S takes equations, but takes advantage of realistic physical approximations to 
1 writ anqwter rw tine and storage requirements associated with the solution 
of the full lavier-Stokes equations. The assumption is wade that a pr inary 
flow direction exists and that diffusion arising fron the rate of change of 
flou in this priw'ry direction can be neglected. The resulting set of equa- 
tions are solved oy an efficient spatial marching pro c ed u re. The code can 
solve the three moment]* equations, continuity equation, energy equation and 
turbulent kinetic energy equation simultaneously. The circ u m f e r e nt ial nonen- 
tiu equation is neglected for tuo-diaensianal cases. For all of the calcula- 
tions in this study, an assumed constant stagnation enthalpy has been substi- 
tuted for the full energy equation, and a simple algebraic turbulence nodel of 
McDonald and Canarata (ref. 9) with Van Driest dapping has been used instead 
of the one equation nodel. 

In each of the test cases the effects of the nuuber and location of 
radial, circumferential and streanwise mesh points has been studied. The 
effects of the nixing length scaling factor delta has also been studied m 
several of the test cases. Because of the many shock-boundary layer inter- 
actions present in inlet flou fields, it is quite difficult to determine the 
proper nixing length scaling. In this study, the scaling was assiwed to be 
the average boundary layer height for the entire surface. On solid surfaces, 
the boundary conditions were zero normal velocity and either zero tangential 
velocity (no-slip) or a slip velocity determined by the law of the wall (wall 
function). On flow surfaces, the boundary conditions were determined by 
extrapolation from interior points along Mach lines which allows flow ami 
shock waves to pass out of the computational domain. A1 ( lengths were non- 
dimensional ized by the cowl lip height. 

Lewis Inlet 40-60 

The first test case considered for PEPSIS verification was the Lewis 40-60 
nixed compression inlet at zero degree angle of attack (ref. 6). The 40-60 
inlet is a Mach 2.5 axisyanetric design with 40 percent of the total compres- 
sion ^occuring upstream of the cowl and 60 percent internally. The inlet has a 
12.5* half angle centerbody, cowl lip aligned with the inlet axis, and inter- 
nal contouring to conduct the flow to the engine face. 

The inlet compresses the flow through a system of shock waves. At the 
design Mach number, the shock generated by the centerbody passes just outside 
the cowl lip. A shock generated by the cowl reflects off the centerbody back 
onto the cowl and again back to the centerbody where the terminal shock is 
formed. To control boundary layer separation, bleed is provided at each of 
these shock reflection locations. Bleed is modeled in the calculations by 
specifying the mass flow through the surface at the appropriate locations; 
uniform bleed was specified, although the code allows for a distribution of 
bleed. 

A two-dimensional calculation was made at design Mach number using an 
89x360 mesh and Reynolds number of 51.0E6. Wall functions were applied to the 
centerbody and cowl surfaces and mesh points were moderately packed near these 
surfaces. A constant step size of 0.02 of the cowl lip height was specified 


2 


and the mixing le ng t h scale delta ws set to 0.03 on the centerbotfy and 0.01 
on the coni. Figure 1 shows Rach wfctr profiles and static pressure distri- 
butions on the cowl and cen t er b ody. The Mach nuhbor profiles indicate that 
PEPSIS goal i tat i rely predicts the proper .-hock formation and reflections. 
Quantitatively, the agre e ment with static pressure is quite good. The upper 
pressure plot shows a steady increase in pres sue on the first part of the 
cowl due to contouring. Ahead of the shock reflection, 3.8 percent of the 
capture flow is bled which low e r s the pressure both analytically and experi- 
nentally. The shock reflected fron the centerbody then causes the pressure to 
rise rapidly. On the centerbody surface, lower curve, the pressure is con- 
stant up to the first shock inpi ngc n e nt where the pressure increases. At this 
shock reflection, 0.28 percent of the capture flow is bled. Downstrean, the 
pressure slowly increases due to the cowl contouring then sharply decreases as 
the flow passes the throat and starts to expand. Just downstream of the 
throat the reflected shock hits and the calculation ends. Figure 2 shows a 
comparison of analytical and experimental total pressure profiles with the 
experiment shown as symbols. The centerbody bleed occurs between the first 
tm> plot locations and the calculation does quite well both fore and aft of 
the :leed zone. Similar plots on the cowl in the high bleed region indicate 
that some refinement of the bleed model is required. 

Figure 3 presents the results of varying some of the input parameters for 
the 40-60 inlet. In figure 3(a) an increase from 40 to 89 radial mesh points 
has provided better detail of the boundary layer structure aft of the first 
bleed zone. A case was also run with the streamwise step reduced from 0.02 to 
0.005 which indicated little change in the profiles throughout the inlet. 
Figure 3(b) shows total pressure profiles for different values of mixing 
length scale delta. The figure indicates that the quantitative results are 
very sensitive to this parameter. An option to automatically determine this 
scaling within the code is being provided. 

Rose Shock-Boundary Layer Interaction 

While the Lewis 40-60 inlet test calculation gave accurate results fore 
and aft of an oblique shock-boundary layer interaction, further detailed veri- 
fication of this interaction was required. The Lewis 40-60 inlet test was 
part of a parametric study of mixed compression inlet performance and data 
from this test was intended only for overall inlet system verification. 
Detailed shock-boundary layer benchmark experimental data was available from a 
study by Rose (ref. 7). His experiment consisted of a cone placed in a blow 
down circular tunnel at Mach 3.88. The shock generated by the cone inter- 
sected the tunnel boundary layer and reflected. The wall region was instru- 
mented with static pressure taps and the boundary layer was measured exten- 
sively through the interaction. 

A two-dimensional calculation was run on an 89x450 mesh with high mesh 
packing on the tunnel wall; there were approximately 30 mesh points in the 
boundary layer. The Mach number was set to 3.88 and the Reynolds number per 
foot was 50.0E6, mixing length scaling delta was 0.10. A boundary layer 
profile for the tunnel wall was required as part of the initial conditions. 
This starting profile was generated by a program from the analysis of Maise 
and McDonald (ref. 10) which is now incorporated in the code as an input 
option. 
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A comparison of the calculated and experimental nail static pressure dis- 
tributions is shown in figure 4. The code correctly predicts the two stage 
pressure rise through the interaction; the first, rather abrupt rise with the 
incoming shock and the second, more gradual rise with the formation of the 
reflected shock wave. Figure 5 shows pitot pressure profiles tArough the 
interaction and a comparison sitb experimental values obtained with a traver- 
sing pitot probe. The incoming shock appears as a discontinuity in the exper- 
imental data, while the shock smearing of the analysis occurs over several 
grid ^ints. Notice that the analysis very well predicts the thinning of the 
bourdary layer through the interaction. 

Rainbird Cone at Angle of Attack 

To test the accuracy of the code for a full three-dimensional flow field, 
the flow about a cone at angle of attack was modeled. The experiment was 
conducted by Rainbird and consisted of a 12.5* half angle cone mounted in a 
wind tunnel which was run at various Mach mmbers, Reynolds numbers and angles 
of attack. A single measuring station was located back along the cone and the 
cone was slowly rotated to obtain circ inferential variations of pressure and 
velocity. Several different Mach nwber and angle of attack combinations were 
tested by the code; this study will present only results of the extreme 
Mach = 4.25, alpha = 15.6* case. While such extreme conditions will probably 
never be encountered by an ax i symmetric mixed compression inlet, the ability 
to accurately calculate this case gives increased confidence in the accuracy 
of the code for much smaller angles of attack. 

The flow in this case figure 6, forms a conical shock of varying strength; 
the windward side being stronger than the leeward. A pressure minimum deve- 
lops about 30* from the leeward ray and an ever thickening boundary layer 
develops from the windward ray to this pressure minimum. Similarly, flow is 
drawn from the leeward ray towards the pressure minimum resulting in a thinner 
boundary layer on the leeward ray. The flow develops a recirculation vortex 
between the leeward ray and pressure minimum as the boundary layers from the 
windward and leeward side interact. 

Several calculations for this case have been run; all of them at a Mach 
number of 4.25, Reynolds number per foot of 51.0E6., delta equal to 0.015, and 
streamwise step size varying from 0.002 to 0.02. Figure 6 shows the results 
of the last PEPSIS calculation for this case. The figure shows Mach number 
contours at a plane cutting through the axis of the cone. The conical shock 
is clearly visible as is the boundary layer structure along the cone surface. 
Qualitatively, the code correctly predicts the shock structure, the thickening 
boundary layer on the cone, the thinning of the boundary layer on the leeward 
side and the recirculation. Figure 7 shows details of this recirculation. 

The upper figure shows secondary velocity vectors and clearly indicates the 
recirculation in the s-shaped profiles near the top. Also evident are 
c-shaped profiles near the side which are also found experimentally. The 
lower figure shows Mach number contours in this region which again indicate a 
recirculation. 

The calculation of figure 7 required an 80x50x220 (radial, circumferen- 
tial, streamwise) mesh. The large number of mesh points were necessary to 
resolve all of the length scales for the important features of this flow 
problem; the creation of the conical shock, the movement of the shock away 
from the cone, the boundary layer and the recirculation region near the cone 
surface. Figure 8 demonstrates the importance of radial and circumferential 
mesh resolution for this problem. The figure shows Mach number contours for 
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the aft quadrant of the flow field for different aeshes at a plane near the 
■easuring station. All of the calculations required 220 streaawise stations, 
with nearly 100 of then located near the cone tip to resolve the formation of 
the shock. Figure 8(a) shows the results of the first calculation with 40 
radial points packed near the cone surface and 19 circurferential points 
evenly spaced at every 10*. The results indicate a thickening of the boundary 
layer on the leeward side, and give no indication of a recirculation. 
Increasing the nuaber of circuaferential p : nts to 37 evenly spaced at 5* 
while aaintaing 40 radial points gives the results of figure 8(b). This cal- 
1 culation gives an indication of recirculation but with little detail. Because 

f the recirculation is caused by interactions in the boundary layer, wore reso- 

* lution was required in the radial direction near the surface. Packing the 40 

radial points closer to the surface, gave the results shown in figure 8(c). 
While resolving details of the surface, resolution of the shock wave was 
sacrificed which resulted in violent post shock pressure oscillations. To 
resolve both the shock and the boundary layer, the radial mesh was increased 
to 80 points. The results shown in figure 8(d) indicate that the calculation 
was then able to resolve boundary layer details without oscillations near the 
shock wave. Figure 8(e) shows one further refinement in which the circum- 
ferential mesh was increased to 50 points which were packed on both the wind- 
ward and the leeward rays. The marked difference in the results of figure 
8(a) and figure 8(e) clearly demonstrates the need to resolve all of the 
length scales present in this problem. 

Figure 9 shows a quantitative comparison of circumferential static 
pressure distribution on the cone surface. The experimental values are given 
by the circles, while the results of a similarity analysis of Jones (ref. 11) 
are given by the squares. The comparison is very good, except near the 
windward ray, where the code underpredicts the pressure rise. There are 
several possible explanations for this discrepancy, including problems in the 
mixing length model near the plane of symmetry, errors in the mixing length 
scaling, or mesh resolution. A comparison of streamwise boundary layer 
profiles indicates that the calculated boundary layer thickness was too small 
for this test case. A thicker boundary layer would give a better agreement 
with experiment because the effective cone angle would be increased. The code 
does correctly model the pressure variations near the leeward ray. 

Lewis 40-60 Inlet at Angle of Attack 

One of the principal uses of the PEPSIS program will be the calculation of 
axisymmetric inlets at angle of attack. A verification test for this problem 
, is presently planned at NASA Lewis using the 40-60 inlet. To demonstrate the 

code's ability to calculate inlets at angle of attack, the 40-60 inlet has 
been calculated at design Mach number at 2* angle of attack and at a Mach 3.0 
, overspeed at 2*. The design Mach number case failed due to subsonic flew on 

the leeward cowl at the first cowl reflection which qualitatively agrees with 
the results of an earlier wind tunnel test. The overspeed case was then cal- 
culated in order to better study the qualitative results of the code, and only 
results of this calculation are reported here. 

For the inlet at design Mach number and at zero degrees angle of attack, 
the shock from the cone passes just off of the cowl lip. For the overspeed 
case, this shock is forced inside the cowl lip. At angle of attack, the shock 
translates forward on the leeward side and falls even farther inside on the 
windward side of the inlet creating a small expansion region on the cowl lip. 
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The calculation used a 49x19x480 mesh and wall function boundary con- 
ditions on solid surfaces. Considering the Rainbird results, boundary layer 
details for this rather coarse aesh are probably not resolved. However, the 
calculation does give the correct qualitative results for the inviscio part of 
the flow field. Figure 10 shows surface static pressure distribution on the 
cowl and cone. The shock waves are clearly shown as dark bands which are in- 
clined relative to the inlet axis. On the cone, the shock systea has moved 
forward on the leeward side and aft on the windward. On the cowl the first 
rather weak shock from the cone tip and the stronger shock from the cone re- 
flection are evident. Also evident on the surface are pressure contours which 
are normal to the inlet axis, particularly near the throat section on the cone 
and along the cowl between shocks. These features are the results of compres- 
sions and expansions present in the inlet due to surface contouring. For 
another perspective, the plots of figure 11 show pressure and Mach number con- 
tours in the flow field at the windward, waterline, and leeward rays. On the 
windward ray, the cone shock falls inside the cowl; while on the leeward ray, 
the shock hits near the cowl lip. Pressure waves from the cowl contouring 
coalesce to strengthen the reflected shock more on the windward than the lee- 
ward. The reflected shock strikes the centerbody further forward on the lee- 
ward than the windward resulting in greater turning of the flow and a stronger 
shock on the leeward. For each shock reflection downstream of the throat, the 
leeward shock is stronger and has moved forward of the windward shock. Also 
evident are the expansions near the throat and its interaction with the shock 
system. Note that the exp’nsion location is fixed geometrically while the 
shocks are not. The Mach number contours indicate this same behavior, but 
also shows the boundary layer variations through the calculation. 

COMPUTING TIME 

Because of the need to perform calculations on very large meshes, a pre- 
mium is put on speed of the computer and the computer code. Table I documents 
the speed characteristics of the PEPSIS code on the NASA Lewis IBM 370 and 
CRAY I. There are two entries for the IBM 370 because the computer can be run 
in two basic modes; either debug, i.e., with an internal symbol dictionary 
(ISD), or in production (without ISD). The ISO imposes a performance penalty 
on the code which is documented in table I. The first two calculations were 
computed in the two-dimensional mode, and its speed advantage over the three- 
dimensional occurs as a result of keeping the whole problem in the core of the 
machine. The figures for the CRAY computer are estimates based on a seven to 
one speed increase over the IBM with ISD. A version of PEPSIS has been run on 
the CRAY to establish this ratio, although the version of the code used for 
this study has not yet been loaded on the CRAY at Lewis. 

CONCLUSIONS 

This paper has presented results from four test cases which are part of 
the verification effort for the three-dimensional supersonic inlet computer 
code PEPSIS. A two-dimensional calculation of the Lewis 40-60 axisymmetric 
mixed compression inlet indicated that the code qualitatively predicts the 
proper shock system formation. Wall static pressure distribution and boundary 
layer profiles fore and aft of a moderate bleed zone compared well with exper- 
imental results. The Rose test case verified the code's ability to quantita- 
tively predict the details of an # oblique shock-boundary layer interaction. 

The Rainbird test case of a 12.5* cone at angle of attack demonstrated the 
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code's ability to accurately model three-dimensional flow fields. The code 
has properly modeled the conical shock wave of varying strength, the cross 
flow boundary layer profiles and the leeward recirculation region. A large 
number of mesh points were required to resolve the details of all of these 
phenomenon; a mesh resolution study indicated that insufficient resolution 
gave either unstable or inaccurate results. Without a highly efficient, high 
speed code like PEPSIS, these detailed calculations s would be impractical. The 
final test case involved the Lewis 40-60 inlet at 2* angle of attack. The 
results correctly model the three-dimensional shift of the shock system within 
axisymnetric inlets at angle of attack. 

The calculation of supersonic inlet flow fields is an extremely complex 
problem because of the many flow phenomenon which occur. A typical axisym- 
metric inlet produces shock waves, expansions, boundary layers and their 
interaction with shock waves, and possible cross flow recirculations. For an 
accurate calculation, the computational mesh must resolve all of the length 
scales for all the phenomenon present in the flow. The Rainbird test case has 
illustrated possible problems associated with insufficient mesh resolution and 
indicated the magnitude of calculations which may be required for proper 
resolution. 

In general, the PEPSIS code gave excellent results for the four test cases 
of this study. Although much work is planned to improve the fluid flow 
physics and decrease the computing time by conversion and optimization on the 
CRAY I computer, it is evident that the PEPSIS computer code is practical for 
three-dimensional supersonic inlet calculations. 
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TABLE I - TEST CASES, MESH AND COMPUTING TIME 


CASE 

MESH 

' 

COMPUTING TIME (MIN) 



IBM 370/3033 

CRAY 1# 



ISD 

NO ISD 


Lewis is 40-60 inlet at 
0° 

89x360 
(32 040) 

5.3 

3.8 

.7 

Rose 

89x450 
(40 050) 

6.6 

4.7 

.9 

Rainbird 

40x19x220 
(167 200) 

42.2 

27.8 

6.0 


80x50x220 
(880 000) 

— 

144.4 

31.7 

Lewis 40-60 inlet at 
0* 

49x19x480 
(446 880) 

117.9 


’6.8 


^ESTIMATED TIME 














Figure 2 . - Centertxxfy boundary layer profiles for Lewis 4D-60 inlet. 
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Figure 3. - Effects of input parameters 
on lewis 40-60 inlet boundary layer 
profiles at rake 2 . 
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Figure 4 - Tunnel wall static pressure distribution for Rose 
shock -boundary layer interaction. 
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Figure 5. - htot pressure profiles for Rose shock-boundary layer interaction. 
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Validation of a Hate Pi-signal Viscous Analysis of 
Aiisjaetric Sopersonic Inlet Flow Fields 
by T. J. Benson and B. H. Anderson 


Rational Aeronautics and Space Aitoinistration 
Lewis Researcb Center 
Cleveland, Ohio 44135 

ABSTRACT 

A three-dimensional viscous narching analysis for supersonic inlets has 
recently been developed. To verify this analysis several bcnchnartr. axisym- 
netric test configurations have been studied and are coppered to experimental 
data. Detailed two-dinensional results for shock-boundary layer interactions 
are presented for flows with and without boundary layer bleed. Three- 
dinensional calculations of a cone at angle of attack and a full inlet at 
angle of attack are also discussed and evaluated. Results of the calculations 
dc n on strate tne code's ability to predict complex flow fields and establishes 
guidelines for future calculations using sinilar codes. 

INTRODUCTION 

Over the past few years, NASA Lewis Research Center has been increasingly 
involved in developing three-dinensional spatial narching analyses for the 
design of supersonic inlets. In support of this work, research was also 
initiated to systematically verify these onputer analyses by a two tier veri- 
fication process. The first tier involves establishing the ability of the 
conputer analyses to predict flow parameters of interest in "generic" con- 
figurations which would be typical of real world designs. A second tier veri- 
fication process was also initiated where the essential elements of the com- 
puter analyses are validated in very simple environments which highlight one 
or more basic or "benchmark" flow mechanisms. This is called "benchmark" 
verification of the computer code. 

A thorough verification process for the PEPSIS computer code is being 
conducted at NASA Lewis. PEPSIS is a three-dimensional viscous marching 
analysis for .upersonic inlets of axisymnetric, rectangular or conformal 
geometry which was jointly developed by Scientific Research Associates and 
NASA Lewis. Details of the algorithmn and code development are given in 
(refs. 1 to 2). Some of the generic verification results have been presented 
in a paper by Anderson and Towne (ref. 3) and at the Lewis Computational Fluid 
Mechanics Workshop in 1980 (ref. 4). Anderson and Benson (ref. 5) present 
some of the benchmark results for rectangular geometry verification, while 
this pao^r presents results for ax i symmetric geometries. Four major test 
case' will be considered in this study; the Lewis 40-60 inlet (ref. 6) at zero 
an$ie of attack, the two-dimensional shock-boundary layer interaction of Rose 
l ref. 7), the cone at angle of attack by Rainbird (ref. 8), and the Lewis 
40-60 inlet at angle of attack. 


AMLYSIS 


PERIS, an acroqn for Parabolic Elliptic Streanwise Inplicit Supersonic, 
is a folly three-dimensional copater design code for sapersonic inlets. The 
analysis retains the g e n eral t hree -dimensional viscous nature of the Ravier- 
Stokes equations, but takes advantage of realistic physical approximations to 
linit computer run tine and storage requirements associated with the solution 
of the full Ravier-Stokes equations. The assumption is nade that a primary 
flow direction exists and that diffusion arising from the rate of change of 
flow in this prin-ry direction can be neglected. The resulting set of equa- 
tions are solved by an efficient spatial marching procedure. The code can 
solve the three nomentia equations, continuity equation, energy equation and 
turbulent kinetic energy equation similtaneously. The circ u m f ere n tial momen- 
ti* equation is neglected for two-dimensional cases. For all of the calcula- 
tions in this study, mi assuned constant stagnation enthalpy has been substi- 
tuted for the full energy equation, and a sinple algebraic turbulence model of 
McDonald and Canarata (ref. 9) with Van Driest dating has been used instead 
of the one equation model. 

In each of the test cases the effects of the n inter and location of 
radial, circumferential and streamwise mesh points has been studied. The 
effects of the mixing length scaling factor delta has also been studied in 
several of the test cases. Because of the many shock-boundary layer inter- 
actions present in inlet flow fields, it is quite difficult to determine the 
proper nixing length scaling. In this study, the scaling was assiMed to be 
the average boundary layer height for the entire surface. On solid surfaces, 
the boundary conditions were zero normal velocity and either zero tangential 
velocity (no-slip) or a slip velocity determined by the law of the wall (wall 
function). On flow surfaces, the boundary conditions were determined by 
extrapolation from interior points along Mach lines which allows flow and 
shock waves to pass out of the computational domain. All lengths were non- 
dimensional i zed by the cowl lip height. 

Lewis Inlet 40-60 

The first test case considered for PEPSIS verification was the Lewis 40-60 
mixed compression inlet at zero degree angle of attack (ref. 6). The 40-60 
inlet is a Mach 2.5 axi symmetric design with 40 percent of the total compres- 
sion # occuring upstream of the cowl and 60 percent internally. The inlet has a 
12.5* half angle centerbody, cowl lip aligned with the inlet axis, and inter- 
nal contouring to conduct the flow to the engine face. 

The inlet congresses the flow through a system of shock waves. At the 
design Mach number, the shock generated by the centerbody passes just outside 
the cowl lip. A shock generated by the cowl reflects off the centerbody back 
onto the cowl and again back to the centerbody where the terminal shock is 
formed. To control boundary layer separation, bleed is provided at each of 
these shock reflection locations. Bleed is modeled in the calculations by 
specifying the mass flow through the surface at the appropriate locations; 
uniform bleed was specified, although the code allows for a distribution of 
bleed. 

A two-dimensional calculation was made at design Mach number using an 
89x360 mesh and Reynolds number of 51.0E6. Wall functions were applied to the 
centerbody and cowl surfaces and mesh points were moderately packed near these 
surfaces. A constant step size of 0.02 of the cowl lip height was specified 
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and tlie nixing length scale delta was set to 0.03 on the centerbody and 0.01 
on the cowl. Figure 1 shows Nach muter profiles and static pressure distri- 
butions on the cowl and centertody. The Nach muter profiles indicate that 
PEPSIS qualitatively predicts the proper .-hock formation and reflections. 
Quantitatively, the ag reemen t with static pressure is quite good. The upper 
pressure plot shows a steady increase in pressure on the first part of the 
cowl due to contouring. Ahead of the shock reflection, 3.8 percent of the 
capture flow is bled which lowers the pressure both analytically and experi- 
mentally. The shock reflected from the centerbody then causes the pressure to 
rise rapidly. On the centerbody surface, lower curve, the pressure is con- 
stant up to the first shock impingement where the pressure increases. At this 
shock reflection, 0.28 percent of the capture flow is bled. Downstream, the 
pressure slowly increases due to the cowl contouring then sharply decreases as 
the flow passes the throat and starts to expand. Just downstream of the 
throat the reflected shock hits and the calculation ends. Figure 2 shows a 
comparison of analytical and experimental total pressure profiles with the 
experiment shown as symbols. The centerbody bleed occurs between the first 
two plot locations and the calculation does quite well both fore and aft of 
the bleed zone. Similar plots on the cowl in the high bleed region indicate 
that some refinement of the bleed model is required. 

Figure 3 presents the results of varying some of the input parameters for 
the 40-60 inlet. In figure 3(a) an increase from 40 to 89 radial mesh points 
has provided better detail of the boundary layer structure aft of the first 
bleed zone. A case was also run with the streamwise step reduced from 0.02 to 
0.005 which indicated little change in the profiles throughout the inlet. 
Figure 3(b) shows total pressure profiles for different values of mixing 
length scale delta. The figure indicates that the quantitative results are 
very sensitive to this parameter. An option to automatically determine this 
scaling within the code is being provided. 

Rose Shock-Boundary Layer Interaction 

While the Lewis 40-60 inlet test calculation gave accurate results fore 
and aft of an oblique shock-boundary layer interaction, further detailed veri- 
fication of this interaction was required. The Lewis 40-60 inlet test was 
part of a parametric study of mixed compression inlet performance and data 
from this test was intended only for overall inlet system verification. 
Detailed shock-boundary layer benchmark experimental data was available from a 
study by Rose (ref. 7). His experiment consisted of a cone placed in a blow 
down circular tunnel at Mach 3.88. The shock generated by the cone inter- 
sected the tunnel boundary layer and reflected. The wall region was instru- 
mented with static pressure taps and the boundary layer was measured exten- 
sively through the interaction. 

A two-dimensional calculation was run on an 89x450 mesh with high mesh 
packing on the tunnel wall; there were approximately 30 mesh points in the 
boundary layer. The Mach number was set to 3.88 and the Reynolds number per 
foot was 50.0E6, mixing length scaling delta was 0.10. A boundary layer 
profile for the tunnel wall was required as part of the initial conditions. 
This starting profile was generated by a program from the analysis of Maise 
and McDonald (ref. 10) which is now incorporated in the code as an input 
option. 
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A comparison of the calculated and experimental wall static pressure dis- 
tributions is shown in figure 4. The code correctly predicts the two stage 
pressure rise through the interaction; the first, rather abrupt rise with the 
incowing shock and the second, wore gradual rise with the formation of the 
reflected shock wave. Figure 5 shows pitot pressure profiles through the 
interaction and a cowparison with experiwental values obtained with a traver- 
sing pitot probe. The incowing shock appears as a discontinuity in the exper- 
iwental data, while the shock swearing of the analysis occurs over several 
grid points. Notice that the analysis very well predicts the thinning of the 
boundary layer through the interaction. 

Rainbird Cone at Angle of Attack 

To test the accuracy of the code for a full three-dimensional flow field, 
the flow about a cone at angle of attack was modeled. The experiment was 
conducted by Rainbird and consisted of a 12.5* half angle cone mounted in a 
wind tunnel which was run at various Mach numbers, Reynolds numbers and angles 
of attack. A single measuring station was located back along the cone and the 
cone was slowly rotated to obtain circtmiferential variations of pressure and 
velocity. Several different Mach number and angle of attack combinations were 
tested by the code; this study will present only results of the extreme 
Mach * 4.25, alpha = 15.6* case. While such extreme conditions will probably 
never be encountered by an ax i symmetric mixed compression inlet, the ability 
to accurately calculate this case gives increased confidence in the accuracy 
of the code for much smaller angles of attack. 

The flow in this case figure 6, forms a conical shock of varying strength; 
the windward side being stronger than the leeward. A pressure minimum deve- 
lops about 30* from the leeward ray and an ever thickening boundary layer 
develops from the windward ray to this pressure minimum. Similarly, flow is 
drawn from the leeward ray towards the pressure minimum resulting in a thinner 
boundary layer on the leeward ray. The flow develops a recirculation vortex 
between the leeward ray and pressure minimum as the boundary layers from the 
windward and leeward side interact. 

Several calculations for this case have been run; all of them at a Mach 
number of 4.25, Reynolds number per foot of 51.0E6., delta equal to 0.015, and 
streamwise step size varying from 0.002 to 0.02. Figure 6 shows the results 
of the last PEPSIS calculation tor this case. The figure shows Mach number 
contours at a plane cutting through the axis of the cone. The conical shock 
is clearly visible as is the boundary layer structure along the cone surface. 
Qualitatively, the code correctly predicts the shock structure, the thickening 
boundary layer on the cone, the thinning of the boundary layer on the leeward 
side and the recirculation. Figure 7 shows details of this recirculation. 

The upper figure shows secondary velocity vectors and clearly indicates the 
recirculation in the s-shaped profiles near the top. Also evident are 
c-shaped profiles near the side which are also found experimentally. The 
lower figure shows Mach number contours in this region which again indicate a 
recirculation. 

The calculation of figure 7 required an 80x50x220 (radial, circumferen- 
tial, streamwise) mesh. The large number of mesh points were necessary to 
resolve all of the length scales for the important features of this flow 
problem; the creation of the conical shock, the movement of the shock away 
from the cone, the boundary layer and the recirculation region near the cone 
surface. Figure 8 demonstrates the importance of radial and circumferential 
mesh resolution for this problem. The figure shows Mach number contours for 
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the aft quadrant of the flow field for different Meshes at a plane near the 
Measuring station. All of the calculations required 220 streawwise stations, 
with nearly 10 of the* located near the cone tip to resolve the formation of 
the shock. Figure 8(a) shows the results of the first calculation with 40 
radial points packed near^the cone surface and 19 circwferential points 
evenly spaced at every 10*. The results indicate a thickening of the boundary 
layer on the leeward side, and give no indication of a recirculation. 
Increasing the nwber of circuwferential points to 37 evenly spaced at 5* 
while Maintaing 40 radial points gives the results of figure 8(b). This cal- 
culation gives an indication of recirculation but with little detail. Because 
the recirculation is caused by interactions in the boundary layer, more reso- 
lution was required in the radial direction near the surface. Packing the 40 
radial points closer to the surface, gave the results shown in figure 8(c). 
While resolving details of the surface, resolution of the shock wave was 
sacrificed which resulted in violent post shock pressure oscillations. To 
resolve both the shock and the boundary layer, the radial mesh was increased 
to 80 points. The results shown in figure 8(d) indicate that the calculation 
was then able to resolve boundary layer details without oscillations near the 
shock wave. Figure 8(e) shows one further refinement in which the circum- 
ferential mesh was increased to 50 points which were packed on both the wind- 
ward and the leeward rays. The marked difference in the results of figure 
8(a) and figure 8(e) clearly demonstrates the need to resolve all of the 
length scales present in this problem. 

Figure 9 shows a quantitative comparison of circumferential static 
pressure distribution on the cone surface. The experimental values are given 
by the circles, while the results of a similarity analysis of Jones (ref. 11) 
are given by the squares. The comparison is very good, except near the 
windward ray, where the code underpredicts the pressure rise. There are 
several possible explanations for this discrepancy, including problems in the 
mixing length model near the plane of symmetry, errors in the mixing length 
scaling, or mesh resolution. A comparison of streamwise boundary layer 
profiles indicates that the calculated boundary layer thickness was too small 
for this test case. A thicker boundary layer would give a better agreement 
with experiment because the effective cone angle would be increased. The code 
does correctly model the pressure variations near the leeward ray. 

Lewis 40-60 Inlet at Angle of Attack 

One of the principal uses of the PEPSIS program will be the calculation of 
axisynmetric inlets at angle of attack. A verification test for this problem 
is presently planned at NASA Lewis using the 40-60 inlet. To demonstrate the 
code's ability to calculate inlets at angle of attack, the 40-60 inlet has 
been calculated at design Mach number at 2* angle of attack and at a Mach 3.0 
overspeed at 2*. The design Mach number case failed due to subsonic flow on 
the leeward cowl at the first cowl reflection which qualitatively agrees with 
the results of an earlier wind tunnel test. The overspeed case was then cal- 
culated in order to better study the qualitative results of the code, and only 
results of this calculation are reported here. 

For the inlet at design Mach number and at zero degrees angle of attack, 
the shock from the cone passes just off of the cowl lip. For the overspeed 
case, this shock is forced inside the cowl lip. At angle of attack, the shock 
translates forward on the leeward side and falls even farther inside on the 
windward side of the inlet creating a small expansion region on the cowl lip. 
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The calculation used a 49x19x480 mesh and wall function boundary con- 
ditions on solid surfaces. Considering the Rainbird results, boundary layer 
details for this rather coarse mesh are probably not resolved. However, the 
calculation does give the correct qualitative results for the inviscid part of 
the flow field. Figure 10 shows surface static pressure distribution on the 
cowl and cone. The shock waves are clearly show.i as dark bands which are in- 
clined relative to the inlet axis. On the cone, the shock system has moved 
forward on the leeward side and aft on the windward. On the cowl the first 
rather weak shock from the cone tip and the stronger shock from the cone re- 
flection are evident. Also evident on the surface are pressure contours which 
are normal to the inlet axis, particularly near the throat section on the cone 
and along the cowl between shocks. These features are the results of compres- 
sions and expansions present in the inlet due to surface contouring. For 
another perspective, the plots of figure 11 show pressure and Mach number con- 
tours in the flow field at the windward, waterline, and leeward rays. On the 
windward ray, the cone shock falls inside the cowl; while on the leeward ray, 
the shock hits near the cowl lip. Pressure waves from the cowl contouring 
coalesce to strengthen the reflected shock more on the windward than the lee- 
ward. The reflected shock strikes the centerbody further forward on the lee- 
ward than the windward resulting in greater turning of the flow and a stronger 
shock on the leeward. For each shock reflection downstream of the throat, the 
leeward shock is stronger and has moved forward of the windward shock. Also 
evident are the expansions near the throat and its interaction with the shock 
system. Note that the expansion location is fixed geometrically while the 
shocks are not. The Mach number contours indicate this same behavior, but 
also shows the boundary layer variations through the calculation. 

COMPUTING TIME 

Because of the need to perform calculations on very large meshes, a pre- 
mium is put on speed of the computer and the computer code. Table I documents 
the speed characteristics of the PEPSIS code on the NASA Lewis IBM 370 and 
CRAY I. There are two entries for the IBM 370 because the computer can be run 
in two basic modes; either debug, i.e., with an internal symbol dictionary 
(ISO), or in production (without ISD). The ISD imposes a performance penalty 
on the code which is documented in table I. The first two calculations were 
computed in the two-dimensional mode, and its speed advantage over the three- 
dimensional occurs as a result of keeping the whole problem in the core of the 
machine. The figures for the CRAY computer are estimates based on a seven to 
one speed increase over the IBM with ISD. A version of PEPSIS has been run on 
the CRAY to establish this ratio, although the version of the code used for 
this study has not yet been loaded on the CRAY at Lewis. 

CONCLUSIONS 

This paper has presented results from four test cases which are part of 
the verification effort for the three-dimensional supersonic inlet computer 
code PEPSIS. A two-dimensional calculation of the Lewis 40-60 axisymmetric 
mixed compression inlet indicated that the code qualitatively predicts the 
proper shock system formation. Wall static pressure distribution and boundary 
layer profiles fore and aft of a moderate bleed zone compared well with exper- 
imental results. The Rose test case verified the code's ability to quantita- 
tively predict the details of an^oblique shock-boundary layer interaction. 

The Rainbird test case of a 12.5* cone at angle of attack demonstrated the 
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code's ability to accurately model three-dimensional flow fields. The code 
has properly modeled the conical shock wave of varying strength, the cross 
flow boundary layer profiles and the leeward recirculation region. A large 
m«ber of mesh points were required to resolve the details of all of these 
phenomenon; a riiesh resolution study indicated that insufficient resolution 
gave either unstable or inaccurate results. Without a highly efficient, high 
speed code like PEPSIS, these detailed calculations^would be impractical. The 
final test case involved the Lewis 40-60 inlet at 2* angle of attack. The 
results correctly model the three-dimensional shift of the shock system within 
ax i symmetric inlets at angle of attack. 

The calculation of supersonic inlet flow fields is an extremely complex 
problem because of the many flow phenomenon which occur. A typical axisym- 
metric inlet produces shock waves, expansions, boundary layers and their 
interaction with shock waves, and possible cross flow recirculations. For an 
accurate calculation, the computational mesh must resolve all of the length 
scales for all the phenomenon present in the flow. The Rainbird test case has 
illustrated possible problems associated with insufficient mesh resolution and 
indicated the magnitude of calculations which may be required for proper 
resolution. 

In general, the PEPSIS code gave excellent results for the four test cases 
of this study. Although much work is planned to improve the fluid flow 
physics and decrease the computing time by conversion and optimization on the 
CRAY I computer, it is evident that the PEPSIS computer code is practical for 
three-dimensional supersonic inlet calculations. 
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TABLE I - TEST CASES, MESH AND COMPUTING TIME 


CASE 

MESH 

COMPUTING TIME (MIN) 



IBM 370/3033 

CRAY 1# 



ISO 

NO ISO 


Lewis is 40-60 inlet at 

0* 

89x360 
(32 040) 

5.3 

3.8 

.7 

Rose 

89x4f^ 
(40 050) 

6.6 

4.7 

.9 

Rainbird 

40x19x220 
(167 200) 

42.2 

27.8 

6.0 


80x50x220 
(880 000) 

— 

144.4 

31.7 

Lewis 40-60 inlet at 
0° 

49x19x480 
(446 880) 

117.9 


16.8 


^ESTIMATED TIME 







Figure 2. - Centerfcody boundary la>er profiles for lewis 40-60 inlet. 
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Figure 3. - Effects of input parameters 
on Lewis 40-60 inlet boundary layer 
profiles at rake 2 . 
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Figure 5. - Pitot pressure profiles for Rose stiock-boundary layer interaction. 





Figure 3. - Effects of mesh resolution on Rainbird con? at angle of attack. 
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